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Abstract -We analyze the stability of micro-domains of ligand-receptor bonds that mediate the 
adhesion of biological model membranes. After evaluating the effects of membrane fluctuations 
on the binding affinity of a single bond, we characterize the organization of bonds within the 
domains by theoretical means. In a large range of parameters, we find the commonly suggested 
dense packing to be separated by a free energy barrier from a regime in which bonds are sparsely 
distributed. If bonds are mobile, a coexistence of the two regimes should emerge, which agrees 
with recent experimental observations. 
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The key step in the recognition process of living cells is 
the establishment of adhesive contacts either between op- 
posing membranes of two cells or between the membrane 
of a cell and the extracellular matrix (ECM). It has been 
shown previously that the organization of bonds within 
domains has strong effects on the adhesion of cells and 
the consequent active response [l]. Most insightful were 
the experiments with cells binding to substrates containing 
ligands organized on a hexagonal lattice of a characteristic 
length between 40 and 150 nm. A length of 58 to 73 nm 
distance between bonds was shown necessary for a suc- 
cessful formation of domains |2], and at distances larger 
then 90 nm, domains would not form 0. 

Instead of using living cells, the so-called bottom up 
approach |4] has been successfully used to elucidate vari- 
ous elements relevant to cell membranes and adhesion [5] . 
The main protagonists of this research are giant unilamel- 
lar vesicles that are functionalized with ligands to interact 
with receptors immobilized on the surface [6]. Depend- 
ing on the density of binders on the substrate and in the 
vesicle, as well as on the intrinsic binding affinity of the 
binding pair, either domains with densely packed bonds 
have been observed to grow radially from a nucleation cen- 
ter, or no specific adhesion was reported. In the context 
of the formation of these densely packed domains, valu- 
able information on their nucleation and the growth [7], 
equilibrium 0, cooperative effects pHTT]. and membrane 
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Fig. 1: The system under investigation: A large patch of a 
bonded membrane that deforms and fluctuates in a harmonic 
potential. Bonds are separated by a distance d. 



roughness [T2HI1] have been discussed over the years from 
both the experimental and the theoretical points of view. 

With the development of experimental methods [15], 
more detailed imaging of the distribution of bonds within 
the adhesion domains became possible. Consequently, 
large domains consisting of sparsely distributed bonds 
have been identified in coexistence with densely packed 
domains p!6llT7] . It was reported that sparse domains may 
become densely packed by a gradual increase of density of 
ligand receptor bonds within an area of the domain of sev- 
eral square microns. However, some sparse domains were 
also found stable on time scales of the experiment (several 
hours) [17 . 

The coexistence between the sparse and dense domains 
driven by membrane- mediated interactions was first shown 
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within an effective model [18]. About the same time, an 
adhesion stabihzed phase separation induced by attrac- 
tive interactions of binders within the same membrane, 
was suggested [19 . A somewhat similar phase diagram 
emerged from considering the interplay of a binding bond 
potential and a non-specific repulsion [20|, but in the ab- 
sence of membrane transmitted correlations. More re- 
cently, a complex phase behavior was suggested for active 
binders and binders of different length coexisting within 
the same membrane [TT1[2T]. Here we show, that the sta- 
bility of sparse domains and their coexistence with dense 
domains emerges from basic principles, at low membrane 
tension. Thereby, the cost of deforming the membrane 
in an effective non-specific potential is balanced with the 
energy gain associated with the formation of bonds. 

The model. — We consider bonds placed on a regular 
square (sq) or central hexagonal (ch) lattice. Thereby, we 
assume that the lateral size of the membrane adhesion 
domain is significantly larger than the size of a lattice 
unit cell. We model the total free energy of a domain 
(containing N bonds) and investigate its dependence on 
the distance d between the bonds. 

The membrane deformation energy (in units of /c^T, 
where /c^ is the Boltzmann constant and T the tempera- 
ture) is described by the Helfrich Hamiltonian p!8 l[2m [22] 
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Thereby, the Monge parametrization is used to represent 
the membrane of a bending stiffness hv as a surface of pro- 
jected area A placed above the substrate at the height 
/i(r), r = {x^y) being the in-plane position vector. The 
first term is the bending energy, whereas the second term 
accounts for the tension a in the membrane. In a simplis- 
tic manner, the last term models the generic membrane- 
substrate interaction potential with the harmonic poten- 
tial of a strength 7 with a minimum at the height ho. In 
the context of mimetic systems, this interaction potential 
encompasses for a number of contributions such as the 
van-der-Waals attraction, or the steric repulsion emerging 
from both repeller and membrane shape fluctuations [22] . 
In the case of cells, numerous other factors associated with 
actin (de)polymerization, active forces, the glycocalix and 
the ECM, may all contribute to this potential depending 
on the cell type and the treatment of the substrate. 

The receptors are modeled as thermalized harmonic 
springs of rest length lo and spring constant A, fixed for 
all bonds on the lattice. When the receptor is relatively 
stiff such as a bulky protein, it is modeled with a very 
large spring constant (A 00). If the receptor is a soft 
polymer, deforming to form a bond, A is set finite. 

Ligands and receptors interact through a square-well 
potential [20 1123] , of a very short range a and depth 65 (1 — 
35 /cbT), the latter associated with the intrinsic binding 



affinity. Thus, the total Hamiltonian H of a domain with 
N bonds situated at position rj^ reads 



n = no^^ - ^0)' + Ne, = NUd + Ne^. (2) 
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Here, /(r;^) is the extension of the j-th spring. When a 
0, obviously h{rl) ^(1*5). Furthermore, Hd denotes the 
deformation energy per bond stored in the membrane and 
all receptors. The last term is the binding enthalpy. 

Free energy. — The stability of the domain is deter- 
mined from the difference AJ^^ between the free energy of 
the domain with N formed bonds J^^, and the free energy 
of the reference state in which receptors and the membrane 
fluctuate freely J^^. Both J^^ and are calculated from 
the partition function comprising all possible conforma- 
tions of N receptors and the membrane. Thereby, the par- 
tition functions of the reference and the bound state {Zub 
and ^5), and hence = \nZub and = \nZb, are as- 
sociated with the conformations in which the membrane 
is outside or within the range a of the square potential at 
the position of all receptors, respectively. The domain is 



stable if AT^ = - J^^ < 0. 

Here we have clearly omitted the change in the mixing 
entropy of binders that could, in principle affect the re- 
sults. However, we assume the adhesion to be mediated 
by large domains of bonds. Such domains typically form 
when the concentration of the mobile ligands is signifi- 
cantly larger than the concentration of the immobilized 
receptors [8^. In the context of the free energy of the 
system, in this regime the mixing entropy will provide a 
constant, proportional to the number of bonds and to the 
chemical potential of the free ligands, but independent on 
the distance between bonds. Consequently, in this regime, 
the mixing entropy will act to re-scale the effective binding 
affinity of a single bond [18] , without affecting the general 
phase behavior of the domains. 

The reference state for N receptors. The partition 
function for the free membrane and N unbound fluctuat- 
ing receptors is, up to the normalization. 
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where C is denoting the result of the functional integral 
V[h'^]. 

In this state, the membrane is free, on average flat, and 
positioned in the minimum of the nonspecific potential 
{hir)) = ho. The fluctuation amplitude is that of a the 
membrane under tension ^24j 
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where the sum runs over all possible wave vectors q with 
q = |q|, set by the system size. The spatial correlation 
function [24l[25] is simply given by 
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COS (q- (rp - r)) 
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where tq and r are arbitrary positions on the membrane. 
In the tensionless limit, Eg = l/(8y^) and G(r) ^ 
— 47r~"'^^^kei(r/^ll), with kei signifying the kelvin function. 

Thereby, = l/(8y^) and = k/j are the vertical 
roughness and the lateral correlation length of a tension- 
less unbound membrane, respectively [22l[2l[[25], setting 
the length and other scales of the system. 

An isolated bond - N = 1. The shape of the mem- 
brane bound to the substrate by one bond is typified by 
the Kelvin function [25]. The associated membrane de- 
formation energy U{o-) = {ho — /i(r5))^/(2I]Q) is quadratic 
with respect to the offset from the minimum of the non- 
specific potential . When A ^ 00 the entire deformation 
is stored in the membrane. If, furthermore, a = 0, one 
finds Uo = {ho — /o)^/2^^, as previously determined [I8] . 

Because both, the membrane and the receptor defor- 
mations are quadratic with the elongation, we map the 
problem of forming one bond to a problem of two one- 
dimensional thermalized springs of stiffness ki and /c2 
(Fig. [2]). The springs are said to interact if their relative 
distance falls within a square-well potential of a (short) 
range a and depth Vo. The the free energy difference be- 
tween the bound and the unbound state 
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is calculated as described previously, by subdividing the 
configurational space of the receptor and the membrane. 
The first term on the right side is identified with the de- 
formation energy of the two springs, characterized by a 
reduced spring constant of the coupled system consisting 
of two springs in series kik2/{ki + /C2) elongated to meet 
the system size L. The second and the third terms are the 
enthalpy gain and the entropy loss due to the formation 
of a bond. Thereby, the last contribution only affects the 
depth of the potential. 

By analogy (Fig. [2]), for a membrane binding to a single 
receptor (indicated by the superscript 1) 
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Thereby, is the total deformation energy associated 
with an isolated bond and sets the energy scale of the 
problem. For a = and A ^ 00, = Furthermore, 

is the effective binding affinity 
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Fig. 2: Mapping of the model for two fluctuating interacting 
springs to the membrane- bond model. 



which is the contribution of a single bond to the free en- 
ergy. In essence, it is the intrinsic binding affinity de- 
creased by the entropic cost related to the change in the 
ffuctuations of the membrane and the receptor upon bind- 
ing. At room temperature and typical parameters this cost 
amounts to several ksT. If mixing entropy of the binders 
is considered in the regime of large domains, the chemical 
potential of free ligands would need to be subtracted on 
the right hand side of the eq. (j8j). 

A domain with N bonds. The free energy of an ar- 
bitrary bond configuration of TV bonds emerges from the 
partition function for the bound system 



N ^ 

Zb = lv[h\r)]l[ I dl{ri)e 
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that accounts for all conformations in which all receptors 
and the membrane are simultaneously within the bond 
potential range a. With the Hamiltonian from eq. ([2]), 
one gets 



2b = 



v/(l + AS]^)^detM 



(10) 



with Mij = {Sij + XG{ri - i^))/(l + AS^) accounting for 
the membrane-coupled deformations of bound receptors 
on positions r and . The free energy difference becomes 



^ 1 
A^^ =n'dJ2 - + 9 1" (det M) . 



(11) 



Similarly to eq. (O, the first term in eq. (pTj) is the 
total deformation energy of the receptors and the mem- 
brane. The second term is proportional to the effective 
binding affinity of a single bond and to the total number 
of bonds within the domain. The last term calculates the 
fluctuation-induced interactions between the bonds, which 
interestingly, can be fully decoupled from other contribu- 
tions. At large bond separations, d Hd ^ and 
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Fig. 3: Free energy density as a function of the distance be- 
tween the bonds for the ch lattice (left), and the respective 
phase diagram (right). All variables are rendered dimension- 
less, a = 0. For a detailed discussion, see the main text. 



Phase diagram. To compare domains with the same 
number of bonds that due to packing, cover different area, 
we analyze the free energy density A/ = AJ^^ /A that if 
negative, signifies stable domains (fig. [3]). Three distinct 
regimes are apparent even though a global minimum is 
found in all cases. This minimum, however, appears at 
such small d that it could be rendered inaccessible by the 
finite size of the binders. 

In the limit when d ^ oo^ the A/ limit is ap- 
proached either from positive values (fig. [3l red-dashed 
lines) or negative values (blue and yellow-dotted lines), 
signifying unstable and stable binding, respectively. In 
the former case (regime {i) in the phase diagram), by a 
maximum at intermediate distances precedes either stable 
(regime (ii)) or unstable domains at large distances be- 
tween the bonds (regime (z)), in which case the domains 
are stable only close to the boundary minimum. The phase 
transition line separating the regimes {i) and {ii) is found 
in the limit of d ^ oo as = 

Following the maximum. A/ curves of the region (ii) 
possess a shallow secondary minimum that is always neg- 
ative, with locally stable domains. The maximum, on the 
other hand, penetrates toward positive values of A/ for 
smaller values of e^, making the domains unstable at inter- 
mediate bond separations. At larger values of e^, the max- 
imum may remain with fully negative free energies, hence 
stable domains can take place for any distances between 
bonds. This is also true for the regime {Hi) in which A/ 
remains negative, despite its monotonous increase. The 
border between the regimes {ii) and {Hi) is given by the 
(dis) appearance of a root in the second derivative of the 
free energy density with respect to the bond distance, and 
has to be determined numerically. 

This phase diagram should be applicable to the situa- 
tion in which the distance between the bonds is predefined, 
such as when one of the binder type is immobilized on the 



substrate [3,6]. The domains should be observable for any 
distance in which the total free energy density is smaller 
than zero. Consequently, the commonly observed densely 
packed domains [16] can be found within our model at low 
free energies. However, when both receptors and ligands 
are mobile, the distance between the bonds within the 
domain becomes a free parameter. In this case, domains 
should be found only at the distances at which the minima 
in the total energy density appear. More specifically, apart 
from the densely packed agglomerates, domains associated 
with the minimum in the region {ii) at intermediate dis- 
tances between the bonds, should be seen. Indeed, coexis- 
tence between densely packed domains and domains with 
a sparse distribution of bonds has been observed recently 
in experiments with mobile ligand-receptor pairs p!7| [26]. 

Optimum deformation. — The understanding of the 
above phase diagram evolves from the analysis of the mean 
membrane shape and the fluctuations amplitude. They 
emerge as moments of the height probability distribution 
p{h{r)) of the membrane within the domain with a fixed 
bond configuration. The latter is a functional integral 
over all appropriately weighed realizations of the mem- 
brane profile 
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Because of the quadratic form of eq. (|2]), p{h{r)) is a 
Gaussian distribution with the expectation value giving 
the equilibrium shape 



{h{r)) =ho- 

and variance S^(r) 
S2(r) - 
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being the fluctuation amplitude. Thereby 
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By setting = and A" = 1 one recovers the results 
presented in previous sections. 

The equilibrium shape can be also determined by direct 
minimization of H from eq. ([2]), with constrained exten- 
sion of bonds, and periodic boundary conditions. Conse- 
quently, the effect of the lattice is explicit, and one obtains 



{h{r)) = ho 



(/> + l/A 



{ho - lo). 



(16) 



Here (j) = a ^ Sq(^^^ + ^q'^ + t) di'^) = 

'^^cos{qr){f<iq^ -\- aq^ + 7)"-*^, with a being the area 



p-4 



Coexistence of dilute and densely packed domains of ligand-receptor bonds in membrane adhesion 












1.002 






1 




(J 


0.998 





sq 








ch 














Fig. 4: For a = 0, stiff bonds, and ^^(cr = 0) = 5.6 feT, 
shapes and the free energy per bond AJ'^/N + Cb for the 
ch and the sq lattice (black and red dotted line, respectively) 
are displayed. The deformation energy per bond Tid on the 
ch lattice is presented with the blue dashed line. Upper 
inset highlights the minimum at = 9-5^||- For tensions 
a = (0.125, 0.25, 0.375, 0.5)^;^ shapes (ch lattice, d = lO^n) 
and the free energy are shown. All curves are scaled by 
that depends on the tension as shown in the lower inset. 



of a unit cell. For the squared lattice, the sums run 
over wave vectors q = (gi,g'2) = {27rzi/d^27rz2/d), with 
^1,^2 being integers, while for the hexagonal lattice, qi = 
27r(2zi + Z2)/{V^d) and q2 = 27rz2/d. 

Combining the total Hamiltonian, eq. (|2]), and the equi- 
librium shape of the membrane, eq. (p!6]) , results in the 
total membrane and spring deformation energy per bond 



ndMr))] 



{ho - hf 
2[^+l/A]' 



(17) 



Thereby, all bonds have the same extension and the mem- 
brane is at h{Y\) = hi), for all r^. This result is consis- 
tent with the first term in eq. (pT]) as well as with the ID 
energy profile found for a membrane deformed by two in- 
finite cylinders |27j, and the potential calculated for the 
interaction between two bonds [18] . 

It is instructive to first analyze the case of stiff bonds 
when A(/i(rb) — /o)^ 0, and H^, eq. ([T7j), depends 
only on the lattice type and d (fig. H]). The optimum de- 
formation energy has a shallow minimum at intermediate 
distances between the bonds, causing the minimum in the 
free energy density. This minimum can be understood by 
analyzing the shape of the membrane. Namely, when the 
bonds are far apart, they act as isolated bonds produc- 








4^11 














sq 


/ / 
f / 








ch 


/ / 

/ / 








1 / / 
f / / 






/ 

/ J 
t // 










4^ 










8^11 



Fig. 5: Maximal fluctuation amplitude for the ch and sq lattice 
(black and red-dashed lines, respectively). Fluctuation maps 
at cr = (horizontal array), and finite tensions at d = 8^|| 
(vertical array) are also shown. Parameters as in Fig. 4. 



ing a local deformation in which the membrane, prior to 
flattening into the minimum of the nonspecific potential 
at /lo, overshoots /iq (fig- SI- As the bonds come closer, 
the overshoots become shared by neighboring bonds, de- 
creasing the overall cost in bending. When the overshoots 
fully overlap, the shallow minimum appears in A/ (shapes 
(c) in fig. |4]). Bringing the bonds even closer, again in- 
creases the bending energy providing an energy barrier 
(shapes (6)). When the overshoots disappear the mem- 
brane starts to flatten between the bonds. Consequently, 
the energy slides towards a boundary minimum 1-id 0, 
as ^ 0. With the increasing tension (at constant 7) the 
overshoots in the shape become less pronounced, but the 
cost for deforming the membrane rise (fig. H]). The sec- 
ondary minimum becomes shallower and appears at larger 
d, shrinking the parameter space for the phase {%%). 

In the limit of d ^ 0, flattening of the membrane caused 
by the large density of bonds makes both the deformation 
energy and fluctuations insensitive to the tension (fig. [5]). 
If, on the other hand, d ^ 00, bonds on any lattice become 
independent from one another and the maximal fluctua- 
tion amplitude of the domain fluctuation profile S(r)^, 
takes the limit Sq. 

Decreasing the spring constant A transmits the deforma- 
tion from the membrane to the receptors, again affecting 
the size of the region (ii) in the phase diagram. For very 
soft bonds (A ^ 0), the entire deformation is stored in the 
springs and (Ji{j^h) ~> /iq)- In this case, the membrane is in 
average flat as evidenced by the mean membrane height h 
within the domain and the bond extensions (fig. [6]). Fur- 
thermore, in this regime, the membrane fluctuates as an 
unbound one, irrespectively of d. 



Conclusions. — We analyzed the properties of large 
adhesion domains forming between a membrane and a flat 
substrate, when the ligand-receptor adhesion competes 
with the nonspecific adhesion. While some aspects of our 
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Fig. 6: Top: Average height of the membrane h (blue) and 
the extension of the bonds (red) as a function of A and the 
corresponding shapes at d = G^y. Bottom: The fluctuation 
amplitude of a bond as a function of A. Full and dotted lines 
are obtained with a = and a = l/{4^']_), respectively. 



model have been investigated previously p!8l[2Q] , our work 
provides a unifying framework within which the stability 
of the domains can be fully explored. Our first conclusion 
is that the energetics of the domains forming on differ- 
ent lattices including the simple hexagonal one (data not 
shown) is qualitatively the same. This suggests that the 
modeling on commonly used sq lattices [24l[28] will well 
reproduce the behavior of domains that most likely form 
on the ch lattice. This result emerges from the decoupling 
of the entropic free energy contributions associated with 
the correlations between bonds and contributions of indi- 
vidual bonds. At room temperatures, the correlation con- 
tributions seem to be small and have no qualitative effects 
on the phase diagram. However, on a level of individual 
bonds, the effect of fluctuations may be significant result- 
ing in considerable differences between the intrinsic and 
the effective binding affinity. Furthermore, another im- 
portant result is the evaluation of several regimes in which 
domains are stable. Consequently, our model provides a 
physical explanation for the recently observed coexistence 
of densely packed and sparse domains of bonds. 

While we have demonstrated the power of our approach 
with a very simple model of receptors and of the bond 
potential, the quadratic nature of the membrane deforma- 
tion energy allows for stability analysis of domains formed 
from bonds with a very wide range of potentials. Such 
extensions of the model may be important for quantita- 
tive comparison with measured data, which is a task that 
may be challenging, both from experimental and theoreti- 
cal points of view. However, in order to fully comprehend 
the stability of adhesion domains, this comparison needs 



to be performed. In this light, the current work could be- 
come the necessary foundation for the understanding of 
the stability of domains in cellular and mimetic systems. 
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